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Abstract 

To Investigate magnetic nozzle processes, a 
theoretical quasl-one dimensional analysis was em- 
ployed. Included in the euialysls are the effects 
of unequal electron and ion temperatures and elec- 
tron thermal conductivity. The acceleration of the 
plasma is controlled by an Imposed magnetic nozzle. 
The parameters that control the flow are the nozzle 
shape, the ratio of kinetic pressure to magnetic 
pressure, and the ratio of the nozzle throat radius 
to the colli slonal mean free path. Results show 
that higher nozzle exit velocities than a corre- 
sponding adiabatic flow are achieved. Also, the 
plasma temperature remains higher throughout the 
nozzle than a corresponding adiabatic flow. 

Introduction 


equations over the control volume shown In Fig. 1 
(this Is a segment of the nozzle). Also, the fol- 
lowing assumptions are made. 

(1) Plasma properties are uniform across the 

nozzle 

(2) No current In x direction, a 0 

(3) Quasi -neutrality, ng = n^ 

( 4 ) Heat flux In r and 6 directions Is 

negligible 

( 5 ) Viscous forces are negligible 

(6) No energy added to flow 

( 7 ) Ratio of electron to Ion mass, m^/m^ = e^, 

negligible ccanpared to 1 

It should be noted that 1 is the basic assumption 
for any quasl-one dimensional analysis. 

Uhder the above ass'jmptions the following 
equations are derived in appendix A. 


Plasma devices such sis the MTO arc(^) and the 
Q machine wind tunnel(^) produce high speed flows 
in a diverging magnetic field. One possible ex- 
planation for the acceleration mechanism in these 
devices is that the diverging magnetic field acts 
as a nozzle. Axial acceleration occurs as a result 
of both the magnetic pressure exerted by the nozzle 
and the conversion of thermal motion into axially 
directed motion. To investigate magnetic nozzle 
processes, a theoretical quasl-one dimensional 
analysis was undertaken. 

The recent analysis of Walker and SeikelC^) 
considers this same problem in a different manner. 
Their analysis considers the flow near the axis of 
the nozzle including the electron thermal conduc- 
tivity. The main assumptions of their analysis are 
that the Ion temperature Is negligible and that the 
Hall parameter (electron cyclotron frequency/ 
electron-ion collision frequency) is a constant, at 
least near the nozzle's axis. 


Vg = = V equal x-canponents of (A4) 

electron and ion velocities 

ApV - (ApV)x=o plasma continuity (AS) 

r, 

AnV - (AnV)x=0 
_ _ ^2 

SV ^ ^ + _i ^ = 0 Plasma momentum (A17) 

dx dx 2piQ dx 


5? 0 

dr y 2 2 p/ dx 
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+ Vei(kTe - ^i) 

V dx dx 2 nV dx V eiv e 1/ 


(A2S) 


The geometry of the problan considered in this 
study is shown in Fig. 1. An imposed axially sym- 
metric magnetic field controls the flow. Any 
change in the magnetic field from induced currents 
is neglected. 

One fairly simple method of aned.yzing this 
plasma flow is by a quasi -one-dimensional analysis 
similar to those applied in ordinary gasdynamlcs. 

In the case of a fully-ionized plasma, however, we 
must consider two species (ions and electrons) ra- 
ther than a single species. Also, since electrons 
are much lighter than ions, the electrons tend to 
have larger random energies and therefore a higher 
temperature than the ions. The small electron mass 
also results in large values of the thermal con- 
ductivity for the electrons. Therefore, in the 
analysis herein unequal electron and ion tempera- 
tures as well as electron thermal conductivity are 
included. 


Quasi -One Dimensional Equations 

The quasi-one dimensional equations are de- 
rived in appendix A by integrating the various 


kT — 

— electron energy (A31) 
A dx 



lx 


+ 


5 dx / 


5 ^ dTg 
2 mg dx 


= 1.866 electron heat flux (A34) 


The bars appearing in the above equations are used 
to denote dimensional variable quantities. The 
subscripts e and i are used to denote electrons 
and ions, respectively. All quantities arejie- 
fined in the usual manner; V is velocity, T is 
temperature, "p is mass^ density, IT is number den- 
sity, p is pressure, B is magnetic field strength 
q is heat flux, and A is the cross sectional 
area of the nozzle at position T[. The quantity 
Vgj^ is the electron-ion collision frequency for 
transfer of mcmentum given by Equa. (A28) in the 
appendix. In appendix A, Equa. (A17) is derived 
without assuming that and ^ are uniform 

across the nozzle. However, in the above equations 
it has been assumed that B^ and Bj. are uniform 
across the nozzle. 


1 


Ill addition to the above equations, we need 
expressions for the total pressure p, and 
For the total pressure we have 

p = rigkTe + nj^kTj = nk(Tj^ + Tg) (l) 


Define the temperature. 




therefore 


P = 


rikT 

o 


( 2 ) 


(3) 


The maenetic field Is assumed to determine the 
nozzle shape. 


dx 


(4) 


where R Is the_^nozzle radius at position x. 

Also, from V . B = 0 and assuming Is uniform 

across the nozzle (see equation inmedlately pre- 
ceeding Eq. (A16)), 

Bj(A = (BxA)x=o (3) 


ii adv I 
5 V dx yjs 


4eH/V= 


dx 


= -1.066 




yq _ q dA 
fcHjjV ’ A dx 


(10a) 


Appearing in Eqs. (7) to (10) are the throat Mach 
number H», the ratio of the throat radius to the 
throat mean free path I»Ael«> magnetic pres- 
sure term G(x), and the square root of the mass 
ratio c. 



We now nondlmenslonallze the equations using 
the following definitions. 



(1-4) 


Q y_Y_x = — T T =^° term Is the ratio of kinetic pressure to 

p*' Y^' 1^' ® ° magnetic pressure at the throat. 



Br 


Br 

^Bx^* 


( 6 ) 


The asterisk subscript is used to denote conditions 
that exist at the nozzle throat. Iftibarred quan- 
tities are dimensionless, T* = (Tq)j^_q and Is 

the throat radius. Also, Vei* hM been defined as 
a function of T* rather than Te* = (Te)x=0* 


If Eqs. (A5) and (3) 
p, FT, and p, Eqs. (A17), 
become, respectively, the 
equations 



are used to eliminate 
(A25), (A31), and (A34) 
following dimensionless 


dT 

^ = G(x) + 

^0 dA 

(7) 

dx 

AV dx 

. il da , 

= 0 
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eM» dx. 
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(The approximation pwaij^ir (eq. (A6)) Is used 
throughout the paper.) Prom Eqs. (7) to (lO) we 
see that the solution will depend on the param- 
eters M)(, (Z«/Xg^ ), and e. Also, we must 

know the area. A, as a*functlon of x. As we 
shall see later, M* will be determined by the 
value of (Tg)xsO' parameter D* determines 

the relative energy density of the magnetic field 
compared to the thermal energy density. Large 
Implies small znagnetlc fields. The parameter 

1 Is a measure of the "cclllsionness" of 
the flow. A large value of this parameter means 
there are many collisions occurring. For the basic 
equations from which the quasi -one dimensional 
equations are derived to be valid, the flow must 
be colllslonal. Therefore, Z«/kej^^ ^ 1. The mass 
ratio c merely depends on the gas considered. 

In Eq. (lOa) the dTg/dx and vq tezms have 
a coefficient proportional to Since this 

quantity is of the order of 10^, the other terms 
in Eq. (lOa) will be small compared to the dTg/dx 
and vq terms. As a result the following approx- 
imate heat flux equation was used. 



If Eqs. ( 7 ) to (IC) are sox ^ for the var- 




lous derivatives the following results are ob- 
tained. 


IX ^ V 
dx 


- — To - Te) + - GV 
Adx\2 ° V 2 


+ a 


dT 


dx 


(16) 


1 ^ 
A dx 


5 * *' o e 


■♦GV 


- l^V^-T 1 + F- M?V^-T 
5 * L3 * 



^ = 
dx 


(17) 

(18) 


^ Isf „ 


( 24 ) 


Equation (23) shows that the critical point occurs 
at a subsonic Mach number. In order for the flow 
to pass smoothly through this singular point, the 
numerators In Eqs. (16), (17), and (19) must go to 
zero when D = 0. The first two terms in the nu- 
merators of Eqs. (16), (17), and (19) are propor- 
tional to dA/dx. Therefore, when dA/dx » 0 
these terms will vanish. The last term In the nu- 
merators of Eqs. (16), (17), and (19) will vanish 
when 0 s 0. As a result, by requiring that da/dx 
=6=0 when D = 0 the solution will be con- 
tinuous at the crltfcdl point. Prom the require- 
ment 9 s 0, we obtain a condition on the heat flux 
at the critical point (see eq. (21)). 


< 1 * 


25 

18.66 




CT, 




(2T, 


“ T(X|j) (25a) 


The reason for the asterisk will be explained 
shortly. 


2^6 


€M» 


fl ^ M^V^T -fOVT -(M^V^-T )o‘ 
Adx**e 


3VT -(M^V^-T )0‘ 
e * o 
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(19) 


Appearing in Eqs. (16) to (19) are the quantities, 

D = - (t„ - I T,jj (20) 


0 = + (To-2Te) 

\\i*j [ 25 Y6T^ ^y2^ 


(21) 


In obtaining these expressions we have used the re- 
sults of Eq. (A28) in the appendix. 


V = — ^ 


n Tn Aei 


Vel, T^/^ ^ ^1* 


(22a) 


Since In Is a slowly varying function we have 

used the approximation. 


AVTg/^ 


(22b) 


The continuity equation has been used to replace n 
In Eq. (22b). 

Critical Point 

The system of Eqs. (l6) to (19) has a critical 
point when D = 0 or when. 




m2 = 




^ kT„ m 2 + ^ c2 

3 o 6 5 


(23) 


where 


For any value of T^^ that gives a physically 
meaningful solution D < 0 Just upstream of the 
critical point (see eq. (20)). Therefore, in order 
for the flow to be accelerating the nvmerator of • 
Eq. (16) must be negative, in magnitude. For the 
values of and l«/Xel^ considered here, the 

predominate term in the numerator of Eq. (16) is 
the first one, which is multiplied by (l/A)(dA/dx). 
Therefore, in order for the numerator to be nega- 
tive, dA/dx < 0. Similarly, Just downstream of the 
critical point D > 0. In this case dA/dx must 
be positive. Since dA/dx goes from negative to 
positive in passing through the critical point we 
have established the critical point as that corre- 
sponding to the minimum area of the nozzle (l.e., 
the throat). It is for this reason that the aster- 
isk subscript is used on the quantities in 
Eq. (25a). Also, since we have nondimensionallzed 
in terms of the throat properties, Tq^ - 1, V* = 1, 
A« = 1 and Eq. (25a) becomes the following. 

25 le , X , » 

V? "iT 

Ifelng Eq. (20) at the critical point we obtain an 
expression for M». 


m2 = 1 - £ T. 
^ 5 ®* 


(26) 


It should be pointed out that the condition 
dA/dx ■ 0 a 0 is not the only one that will make 
the numerators of Eqs. (16), (17), and (19) vanish. 
However, based on ordinary gasdynamlc flows, it 
appears physically meanln^^ul that the critical 
point should occur at the throat. Also, the re- 
quirement 0 ■ 0 establishes a boundary condition 


on the heat flux in terms of T, 


e»* 


Hozzle Area Ratio 


So far it has been established that the equa- 
tions describing the flow have a critical point. 
Further, this criticsil point has been fixed at the 
nozzle throat and yielded a boundary condition on 


3 



the heat flux. This information was obtained with- 
out specifying the cross-sectional area variation 
of the nozzle. To integrate the equations, how- 
ever, we must si5ecify how the nozzle area varies 
with X. 

The nozzle radius dependence on x was 
assumed to be the following. 

R a aj?* + (27) 



where Z* is the throat radius and a and p are 
constants. Using Eq. (.77) the area can be calcu- 
lated. In dimensionless form it is the following. 


A = (AgxP + 1)‘ 


(28) 


where. 


Ag = aZ* 


P-1 


From Eq. (20) the following is obtained. 


A dx 




Also, from Eqs. (4), (s), and (27) to (29), 

,P-1 

~? 




(AgxP + !)■ 


(23) 


(30) 


(31) 


and using Eq. (13), 


G = 


2p(p_^:r^ 1 

p* [p - 1 ® 


,2p-3 


(AgXP + 1)' 


(32) 


Dependence of Solution on Electron Temperature 

Considering Eqs. (16) to (19) and (25) to 
(32), we see that the solution will depend on the 
following quantities; c, (Z*/\el*)> Ag, p, and 
Tg . The first three of these parameters have al- 
reSdy been discussed. The quantities Ag and p 
merely determine how fast the area changes with x. 

Since is a boiuidary condition on the 

electron temperature it might at first seem tliat 
solutions should be obtainable for all values of 
Tg^ from 0 to 1. However, it was found that for 
given values of the other parameters a physically 
meaningful solution was obtained for only one valua 
of Tg^. By a physically meaningful sol’jtion is 
meant one in which Tq ^ Tg for large velues of x. 

Asyir.t’totlc Solution 

For very large x the expected results are 
Tq# Tg, q -* 0 and V — V^, where V is soro* 
asymptotic limit. However, Eqs. (16) to (19) will 
not produce such results. This can be seen by con- 
sidering Eq. (19) for large x. In this case 


where Eq. (18) has been used. From this expression 
we see that q will be Increasing while Tg is 
decreasing for large x. Such a result is not 
physically meaningful since both Tg and q 
nhould approach zero as x approaches infinity, 
'.’herefore, the integration of Eqs. (16) to (19) was 
■:ut off when q started to Increase. This oc- 
curred for values of the area ratio, A of the 
order of 100. A discussion of how tills result in- 
fluences the conclueic>us of the analysis is pre- 
seupeu lu the RESULTS section. 

Adiabatic and Isentroplc Solutions 


It is of Interest to compare tlis solution o* 
Eqs. (16) to (19) with the adiabatic and isentroplc 
flows for a similar nozzle. To obtain the adia- 
batic solution we merely let q a 0 in Eqs. (7) to 
(10). In this case Eq. (lO) is not used and 
Eq. (9) (electron energy) - is decoupled from 
Eqs. (7) and (8). Equation (8) can be Integrated 
to obtain 


To 


— (1 - V^) + 1 

3 


(33a) 


We also obtain from 
results for dV/dx 


dx 


Eqs. (7) to (9) the following 
and dTg/dx; 


. V 


|rjV + -2. 

A dx 


■ Tc 


(34a) 




GV + i — 

A dx 

Mi'V^ - Tg 


(35a) 


Tills system of equations has a critical point vtnen, 


- Tg =. 0 


or when 



Thi.s is the ordinary gasdynamic situation vliere the 
critical point cccurs at M * 1. Similar to the 
case for q / 0, the critical point must occur at 
the nozzle throat in order for the flow to pass 
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amo'thly through this point. Therefore, since M* 
is tho throat Mach number, = 1 and Eqs. (35a) 
to (.’ha) become, respectively 


To = 7 (1 - V") + 1 



(33b) 


(34b) 



GV + - — 

A dx 

2 

V -T„ 


(35b) 


Equation (33b) yields 
ratio that can be attained 
a monatomic gas. To reach 
the flow is expanded until 
from Eq. (33b) 


the maximum velocity 
in an adiabatic flow for 
the maximum velocity, 

Tq -♦ 0. Therefore, 


(VraEoc)ad ” ^ 


(37) 


For isentropic flow, q = 0, and no current may 
flow in the plasma^ Tlje magnetic field term G 
results from the J x B term in the plasma momen- 
tum equation. Therefore, since J -» 0, G -• 0 as 
well. Under these conditions the usual one- 
dimensional isentropic flow equations for a mona- 
tomic gas (t = 5 / 3 )) are obtained. (4) The temper- 
ature and velocity are related by Eq. (33b) and the 
flow Mach number is given as a function of the area 
ratio A by the following expression. 

W* + 6M^ - 16 AM + 9 = 0 (38) 

Since M* = 1> 

V = (39) 

Using Eq. (39) in (33b) we obtain Tq as a func- 
tion of M. 

Tc. = — (40) 

3 + 

Equation (33) can be solved for M(A) and then 
used in Eqs. (39) and (40) to obtain V and Tg 
as functions of the area ratio A. 


points very close to x » 0 did not greatly effect 
the solution obtained. For all the cases pre- 
sented here the following starting value were 
used; x^ = 0.001, V* = 1.0001, and Tg^ -1.0. An 
iterative process was necessary to establish the 
proper value of the initial electron temperature 
Tg*. The critical heat flux q*, and Mach number 
were csilculsted using Eqs. (25b) and (26), re- 
spectively. Solutions were obtained for various 
values of the parameters 3# and (l»/Xei*) In 
argon (c = 3.71x10"^). 

In Fig. 2 the complete solution is compared 
to the adiabatic and isentropic solutions for 
Ag = 0.02, p = 2, P* = 5, and (l*Ael#) “ 5. The 
velocity and temperature profiles for the adiabatic 
and isentropic cases are nearly the same. That is 
why they are shown together in Fig. 2. There are 
two significant things to notice from these re- 
sults. First of all, velocities higher than the 
adiabatic limit can be attained when q / 0. Sec- 
ondly, the electron temperature, and likewise the 
total temperature, remains relatively high through- 
out the nozzle. The cause of this temperature re- 
sult is the large thermal conductivity for the 
electrons . 

Changing P# while keeping the other param- 
eters constant results in only a slight change In 
the velocity and heat flux profiles. Increasing 
P« produces a smeLll increase in the velocity. 

The temperature profiles show a greater change with 
P«. Figvure 3 Illustrates this effect. Larger 
vsdues of P^ result in higher temperatures 
throughout the nozzle. Such a result is expected 
since large values of P« mean that the thermal 
energy density is large compared to the magnetic 
energy density. 

The effect of the collision parameter !«/ 

Xei , is more pronounced than P«. Increasing the 
value of l*Aei* produces smaller values for V, 
Tg, Tg> and q. Figure 4 shows this effect on V, 
and Tq. Large values of 1*/Xei« imply a more 
collislonal flow. As the flow becomes more colli- 
sional the electron thermal conductivity is de- 
creased. Therefore, the flow approaches the adia- 
batic limit and the resiilting lower values for V, 
Tg, Te, and q. 

In Ref. 5 the plasma properties were measured 
in the magnetic nozzle of a low power MTO arc. It 
is difficult to compare these experloenteLl results 
with the theory since the values of P^ and l«/ 
^ei* appropriate to the experiment are not known. 
However, there is qualitative agreement between 
the theory and experiment for the electron tempera- 
ture. The experimental electron temperature de- 
creases much slower than that in an adiabatic ex- 
pansion which is in agreement with the theory. 


Results for the adiabatic (eqs. (33b) to (35b)) 
and Isentropic (eqs. (38) to (40)) cases will be 
com[iared to the complete solution in the following 
section. 

Results 

Equations (IG) to (19) were integrated using 
the Runge-Kutta method. The integration was 
started near the throat of the nozzle ( x = 0) . 

Since the throat is at the critical point it is not 
possible to begin the integration at exactly x = 

0. However, it was found that different starting 


It has already been pointed out that the as- 
ymptotic behavior of the heat flux equation is not 
ph^icaliy meaningful. This raises two questions. 
Is the heat flux equation correct? If it is in- 
correct, how are the results effected. The first 
question can not be answered without experimentally 
testing the electron heat flux equation used in 
this analysis (eq. (A34)). However, a qualitative 
answer can be given for the second question. As 
long as electron beat flux, q, exists, no matter 
what equation we use to describe it, the conclu- 
sions about higher temperatures and velocities 


5 


than a corroBfxjndlng adiabatic flow (q = O) will 
still be valid. This is so because the presence of 
heat flux gives the electrons another mechanism for 
transporting their energy down the nozzle and 
therefore producing higher temperatures. This addi- 
tional transport of thermal energy down the nozzle 
eventually ends up being converted into directed 
energy (vVS). As a result higher velocities than 
the adiabatic case (q = O) will always be attained. 

Conclusion 

Two significant results were obtained from the 
quasl-ono dimensional analysis of a magnetic nozzle. 
First of all, higher velocities than the adiabatic 
limit can be attained in the nozzle. Secondly, the 
electron temperature and therefore the total tem- 
perature remains high throughout the nozzle. Both 
of these results occur because of the Inclusion of 
the electron thermal conductivity. The experimen- 
tal electron temperature measurements of Ref. 5 in- 
dicate that the temperature remains high throughout 
the nozzle. This result agrees qualitatively with 
the theoretical result. 

Appendix A 


Assuming Jv = 0 and charge neutrality (nj^ - 
Ug = n), Eq. (A3) yields, 

Ve = Vi = V (A4) 

where V is the x component of the plasma mean 
velocity. Also, Eq. (A2) for the ions and elec- 
trons together with the result (eq. (A4)) can be 
used to obtain the continuity equation for the 
plasma. 

ApV - (ApV)x=o (A5) 

where 

P = Pg + Pi = ®i*^l + ®e”e “ "^®1 ®e^ 

and since m^ « mi, 

p wpi = min (Afi) 

The steady-state plasma momentum equation is, 
(in MCS units) neglecting the gravitational term, 
the following (eq. 17.14 of ref. 6). 


Quasi -One Dimensional Equation s 

The steady-state quasi-cne dimensionsd. equa- 
tions are derived by integrating the plasma equa- 
tions over a control volume (which is a segment of 
the nozzle) such as that shown in Fig. 1. For ex- 
ample, consider the steady-state continuity equa- 
tion for the ions or electrons. 

— (nsUsi) = 0 (Al) 

OXi 


pUi 




dxi dxj 


- JlBmClmJ - oEj 


0 


(A7) 


Vfhere, Ui is the plasma velocity, Pij is the 
stress tensor, p is the scalar pressure, Jj is 
the total current density, is the magnetic 
field strength, cimj is the alternating tensor (or 
Levi-cevlta denslt^, c is the net charge den- 
sity, and E. is the electric field. Since n m 
ni, 0=0. Also, neglecting the stress tensor the 
following is obtained from Eq. (A7). 


The subscript s denotes either electrons or ions. 
Also, cartesian tensor notation is used. 

Integrating Eq. (Al) over the control volume 
yields 

^ (nsUsi)dr = 0 

OX 

Now apply the divergence theorem and assume nsUsj^ 
is uniform across the nozzle. 

“ AngVg - (AngVg)^_Q = 0 (A<^) 

In obtaining Eq. (A?) the condition Ugj^hi = 0 on 
the nozzle's lateral surfEice, was used, where h 
is the unit vector perpendicular to the control 
volume surface 3. The quantity Vg is the x 
component of the velocity of species s(Vg > O). 

Equation {k?) holds for both the Ions and 
electrons. Using Eq. (A") for the ions and elec- 
trons the following expression for the current den- 
sity In the X direction is obtained. 

J., = e(niVi - neV^) = i (AniVi - AneVg)^^Q(A3) 
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pUi 


dx^ Sxj 


Jl®mcimj 


(A8) 


Using the steady state plasma continuity equation 
(eq. 17.9 of ref. 6) 


a(pUt) 

dx^ 


0 


(A9) 


in the first term of Eq. (Ab) and Meucwell's equa- 
tion. 


aiBj 

Hii ■ >oh 


(AlO) 


for J^, the following result is obtained. 


^ foUiUj - ^ [p + ^ 


0 (All) 


where the fact that dB^^/dx^ ■ 0 has been used. 


If Eq. (All) is now Integrated over the con- 
txx)l volume shown in Fig. 1 and the divergence 
theorem applied, the following is obtained. 



where e Is the magnitude of the electron charge. 


since the nozzle surface S' Is assumed to be de- 
termined by the magnetic field, B^hi = 0 and 
Uihi =0 on S'. Therefore, 


Since B^h| > 0 on the nozzle surface this becomes 


pVUj - ^ t 
4o 


•Ao=A(x=0) 








The quantity 5xj = 0 for x J and 5xJ “ ^ 
if X - J. 

h dS' = - ^ dx' (A13) 

^ dx' 

Using Eq. (A13) and the assumption that the plasma 
properties are uniform across the nozzle gives the 
following for the x component of Eq. (A12) 

ApV^ + pA - I^ApV^ + - J 


^ (B? - 


Bi)dA' 


■ iT ‘ 

*^o •'q 

where the fact that B„ = 0 tias been used. 

e 

Now take the x derivative, using Leibnitz's 
rule to differentiate the Integral terms. 

pV ^ ^ ^ ^ dA' 

dx dx 2upA 5x 

- dA' + 2B^| ^ * 0 (A14] 

2Uo^ Jq hx |r=R 3x 

where the result ApV = constant (eq. (.\S)) has 
been used. The last term in this expression cm 
be rewritten using Maxwell's equation 

.. SB. 

V • B = — i- ■» 0 (A15] 

5xi 


Integrating Eq. (AIS) over the volume using the di- 
vergence theorem yields the following 

f df ■ B^b^dS ■ 0 




B^dA; - 0 


Taking the x derivative yields 
« I fiA /A(x) 


I dA . . ^x 

dx dx 


Substituting Eq. (A16) in (A14) results in the 
plasma momentum equation. 


dx dx 2i^ Uj Sx 


0 (A17) 


It should be noted that the plasma properties have 
been assumed constant across the nozzle in obtadn- 
ing Eq. (A17). However, the magnetic field has not 
been assumed uniform across the nozzle. If we 
assume that B„ is uniform so that Bx “ Bx|j.^,- 
then the B^ term in Bq. (A17) vanishes. 


The steady-state plasma energy equation is the 
following (eq. 17.20, rtf. 6) 


L + -P — +P«, — i + — i 

dxj^ 2 dxj dxj^ dxj 


Jl(Ei + UiB„cim) (A18) 


Where q^^ is the toted heat flux defined with re- 
spect to the plasma flow velocity. If the mcnentum 
equation (eq. (A7)) Is multiplied by Uj and the 
result added to Bq. (AlC), the following Is ob- 
tained. 


IT * t ^ 


+ ^ (A19) 

dXi 

Now use the steady-state continuity equation, 

d(pUi) ^ ^ 

1 _ A f ^ on\ 


to rewrite the first term in Eq. (A19), 

i!;’ [““‘{f * t o)] * ^ ^ 
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Using Maxwell's equation 


* 

(VxB)i - WoJi (A25?) 

OX^ 

the term can he written aa 



To complete the set of quasi -one dimensional 
equations we need the steady state electron energy 
equation and the heat flux equation. The heat flux 
q that appears in Eq. (A25) is the sum of the 
electron and ion heat fluxes, however, since 
q^ is of the order (ref. 6, p. 62, 

p. 155), the ion heat flux q^oni vill he neglected. 
Therefore, <lx “ *le * 

The s.eady-state electron energy equation is 
the following (ref. 6, eq. (12. 4a)) 


since > 0 In steady state conditions. There- 

fore, Eq. (A21) becomes the following. 


dx^ I y2 2 p/J dxj Sxj 

. 


€,i„ - 0 (A23) 


If we now neglect the viscous term and the Poynting 
vector term E^B„ and integrate the result over 
the control volume in Pig. 1 using the divergence 
theorem, the following is obtained. 


u, ^ /i p.\ . 6 p, ^ . p, 5!li . 

^ Sxj^ \2 / 2 dxj^ &xj^ dxj 


3m^m 


'e“ion 


t“e ♦ “ion)‘ 


“eVei^CTi - T.) (A26) 


The X component of the ste.. .-state electron heat 
flux equation is also obtained frosi Ref. 6 
(eq. 12.9) 


d, 


•tex 7 ^ 7 ^ ^^x ^ 2 „ 

* 5 Sxj 5 Sx^ ” ^ 


1 hx 



0 (A24) 


Now assume the plasma properties are uniform across 
the nozzle and that qj. and q^ are negligible. 
Then, using the fact that Vjh^ -0 on S' 

Eq. (A24) becomes 




where Ug and U_ have been usumed small com- 
pared to Ujj ■ V ti.e., U^ m v"). How take the x 
derivative and <jse the continuity equation 
(eq. (A5)). 





^ i «eiVlax - ^ei^y 


(A27) 


In obtaining the right hand sides of Bqs. (A26) sad 
(A27) results of Ref. 6 given on pages 63 and 150- 
152 were used (the diffusion velocities sre assumed 
to be small). The quantity is the electron- 

ion collision frequency for transfer or momentun 
given In MC3 units by the following expression 
(ref. 7, p. 252), 


In A«4 

Vei • 3.62x10 ** n^ 


3/2 

m • 


A.1 ■ 1-25x10’^ 


(A28) 


(A29) 


It should be noted that the heat flux vector used 
in Eq. (A25) is defined in terns of the mean gas 
velocity while the heat flux vector used in 
Eqs. (A26) and (A27) is defined in terms of the 
species mean velocity u«^. However, in tue prob- 
lem being considered u^^ ■ uionx “ Ux so that 
the qx of Eq. (A2S) is equivalent to qex of 
Eqs. (A26) and (A27). 


pV 



(A25) 


Equation (A25) is the steady-state plasma energy 
equation. In obtaining Eq. (A25) the electromag- 
netic energy transfer given by the Psynting vector 
has been ignored. 


Now convert Eq. (A26) into the quasi -one di- 
mensional fora. First of all, neglect the viscous 
term, Peij. Also since (me/mion)« ■e*lon/ 

“e “lonr • "eAlon* Therefore, Bq. (A26) can be 
written as the following. 
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Whore, - (n»c/™lon)' Integrating fiq. (ASO) over 
the contlMl Volume In Pig. 1 and using the diver- 
gence theorem results in the following. 


IS ii . 5- /i kT„\ . i- iii 

V dx dx \Z / nV dx 

- - ^ Ve ,k(Te - T^) - -1 ^ (A3l) 
V i ^ A dx 

Equation (A3l) is the quasi -one dimensional electron 
energy equation. 



[i \T [t * 


♦ (f Pe«ei ♦ ‘leij 


.IhidS' 




3c^neVeik(Te 



Now consider the quasi-one dimensional fora of 
the electron heat flux equation. Neglect the stress 
tensor terms in Eq. (AS?) and rewrite as follows. 






♦ 


7 ^ 

5 **®1 &xi 


~ ‘lei®m*lmx 
“e 


■ -1.866 v« q, 
®l ^x 


Using the relations dT A dx’, h^^dS' « -(dA/ 
dx’)dx', Ug h^ » 0 on the surface S’ and assum- 
ing that qji anu q^ are negligible the following 
is obtained 


f. 


A(x) 



Tt 


PeV 


♦ qxl 





3c^neVejk(Te - Ti)Ia dx' 


Now use the assumption that the plasma properties 
are uniform across the nozzle. 


Integrating over the control voIubm in Pig. 1 and 
applying the divergence theorem yle\ds the follow- 
ing. 


•X) 0 

‘ 4 * I / ^ * ’«i 



(A32) 


{l * ■‘-j ■ *"f ["''' * "‘I ■ / ^ 


dx* 


Take the x component of Eq. (A3S). Using the 
approxlawtlons that qgj. and qe^ are negligible 
and that the plassu properties are uniform across 
the nozzle yields 


“ ^ VV - (lxAV)x.o *r jT ^ ^ 


where ue^ and Ue^ have beet, neglected in com- 
parlsun to u«^ » V. Taking the x derivative and 
using Pg * ne*^e yields the follow- 

(A33) 




■/ 


-1.866 1 VeiqxA dx* 
^0 


n..V kT \ ♦ - VkTg file - -’^\vg,k(T^-T. ) 

dxy" J dx dx 

Using the electron continuity equation to replace 
dn,,/dx results In the following. 


where the subscript "e” on q^ has been dropped 
since qg^ m q^. In obtaining Eq. (A33) the re- 
sults dr ■ A dx* and Ue^h^ *0 on 3* have 
been used. Taking the x derivative results In 
the quasi -one dimensional electron beat flux equa- 
tion. 


3 



V +. q r J l 4- k 

dx ^ [^A dx 5 dxj 2 ny, ux 

- -1.B66 vel^x (A3*) 

The electron pressure, p . has been replaced by 
nokTg « nXTe In Eq. (A32). 

The complete set of quasl-one dimensional 
equations for p, V, Tj, T , and q^ are given by 
Eqs. (A5), (A17), (A25), (A31), and (A34). 
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Figure 1. - Magnetic nozzle. 
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(A) VaOCITY PROFILE. 

Figure 2. - Comparison of solution with adiabatic and 
isentropic solutions. Magnetic nozzle parameters, 
p - 2 and A^ - 0.02, ratio kinetic to magnetic pres- 
sure Oo ' 3. collision parameter Ul Xg| > 5. 
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(B) TEMPERATURE PROFILES. 


Figure 2. - Continued. 
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(C) HEAT aUX PROFILE. 
Figure 2. - Concluded. 
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Figure 3. - Effect of kinetic to magnetic pressure ratio, 
Oo, on total temperature. Magnetic nozzle parameter 
p-2 and A 5 -O. 02 , collision parameter l«/Xgj -5. 
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(A) VELOCITY PROFIlf. 

Figure 4. - Effect of collision parameter on velocity and 
total temperature. Magnetic nozzle parameter p - 2 
and Aj - 0 . 02 , ratio of kinetic to magnetic pressure 
Do " 5. 
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(B) TEMPERATURE PROFILE. 


Figure 4. - Concluded. 
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